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Abstract. A recently proposed step-by-step procedure, to merge the low-energy physics of the 
TT-bonds electrons of graphene, and quantum field theory on curved spacetimes, is recalled. The 
last step there is the proposal of an experiment to test a Hawking-Unruh effect, emerging from 
O ' the model, that manifests itself as an exact (within the model) prediction for the electronic 

"sJ ', local density of states, in the ideal case of the graphene membrane shaped as a Beltrami 

5_j ■ pseudosphere. A discussion about one particular attempt to experimentally test the model 

&J0' on molecular graphene is presented, and it is taken as an excuse to solve some basic issues 

that will help future experiments. In particular, it is stated that the effect should be visible on 
generic surfaces of constant negative Gaussian curvature, that are infinite in number. 

CN ■ 1. Introduction 

'sj" , Recently (see [1], [2] and, for a review, [3]), in order to realize in practise a quantum Dirac field 

^^ [ theory on curved spacetimes, there has been interest in exploiting certain special features of the 

low-energy physics of the vr-bonds electrons of graphene, Weyl symmetry [4J being among them 
[1]. The clearest signature of quantum field theory (QFT) on curved spacetime is the Hawking 
effect [5]. In [2] the best conditions for the simplest realization of the latter on graphene were 

p\ . proposed. 

^ I On one side, we have graphene, a one-atom-thick membrane of carbon atoms arranged in a 

hexagonal lattice, see Fig. [H first discussed theoretically in [6l [7], then found experimentally 
in [^. Its very peculiar transport properties are summarized by the following low-energy tight- 
binding Hamiltonian, see, e.g., [9], valid near the Fermi points k^ = (±—^,0) in the Brillouin 
zone, see Figs. [2]and[3l 

O-f-f p 

p 

where t ~ 2.8eV, is the electron energy for hopping to the nearest-neighbor, i ~ 1.4A is the 
carbon to carbon distance in the hexagonal honeycomb lattice, vp = 3ti/{2h) ~ c/300 ~ lO^m/s, 
and with a = ((Ti,o"2), a* = (— (Ti,(T2), <7j the Pauli matrices and ip± = {a±,b±) Dirac spinors, 
where a± = a{k^) {b± = b{k^)) is the annihilation operator for the electron in the La {Lb) 
sublattice, evaluated at the Fermi points. 
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Figure 1. The honeycomb graphene lattice. The basis vectors, oi = ^{V^^'i), ^2 = f (\/3; —3), 
hnk sites of the sublattice La, but cannot reach the sublattice Lb- For the latter sites the extra 
vectors, si = £{0, —1), S2 = n (v^j 1); ^3 = 2(~^' ■'■)' ^^^ necessary. 




Figure 2. The dispersion relations of 
graphene, E{kx, ky). The energy E is in 
units of the hopping, t ~ 2.8 eV. The 
conductivity band (red) and the valence 
band (green) touch in six points (where 
E = 0), reflecting, in the Brillouin zone, 
the hexagonal symmetry of the lattice. 
Only the two points kP are inequivalent. 




Figure 3. The linear dispersion relations near one 
of the Dirac points, showing the typical behavior 
of a relativistic-like system (the "vi^'-light-cone" ) . 
Note that the energy range is of the order of 
one unit of the hopping energy t, i.e., here E £ 
[— 3eV, +3eV]. Within this range, the continuum 
field approximation is tenable. 



On the other side we have QFT in curved spacetime. This is by now an old theory [5], 
with no ultimate experimental evidence, hence with many open problems, and interesting 
issues to investigate in relation to the practical implementation of a fundamental theory. The 
open questions range from the structure of the quantum vacuum, and the related role of the 
inequivalent quantizations (think, for instance, to the querelle between the young and the old 
Hawking, about the information loss in the presence of a black hole), to the very meaning 
of particles in a curved spacetime. There are then open questions relative to the Hawking 
phenomenon itself, predicted but to-date impossible to detect from astrophysical sources, and the 
associated black-hole thermodynamics. These questions, although old, are all still on the table 
and permeate the contemporary research on high energy theory, as, for instance, one direction 
pursued these days is to see gravity, and the spacetime itself, as an emergent phenomenon. 

Graphene appears to be a system where many of the issues above could be experimentally 



addressed. In the following, I will propose arguments and experiments in favor of this. 

2. The step-by-step merging 

The preliminary step to start the merging is to consider the dimensionality of the problem, that 
is (2+1) dimensions. I will take this duly into account. For instance, it is by now well known 
that many features of (3+1) dimensions need to be changed on the gravity side [TO], [11], [12], 
|13] . that a black-hole is possible p3], and that certain features, on the matter side, such as the 
"statistical swapping" [15j, can take place. 

• The 1*^* step is to take time on board, hence we move from the Hamiltonian ([1]) to the 
corresponding action 

A = ifiVF I d^x'4^ -i^da V' , (2) 



where a = 0,1,2, 7*^ = 173, 7^ = ia2, 7^ = — «o"i, and, since we do not consider phenomena 
mixing the two Fermi points, we focus on one point only, e.g., V = "0+- When the graphene 
sheet is curved, and within certain limits f3|, the action to consider is 

A = ihvF J d?x^ V57^V/,V , (3) 

where V^ = 9^ + Jl^, and Q.^ = ^uj/^'^Jab is the gauge field able to take into account intrinsic 
curvature. Here a; " is the spin connection, and Jab are the generators of local Lorentz 
transformations, with fi, ... responding to diffeomorphisms, and a,b, ... flat indices. Within the 
limits indicated in |3] and [16] (essentially, in the absence of torsion and of strong scatters J17j , 
and for small curvatures, i.e. large radius of curvature w.r.t. the lattice, r > i), this action well 
describes the physics of the very long wavelength it electrons on curved graphene sheets. The 
only effects it needs to take care of are the ones induced by intrinsic curvature, i.e., inelastic 
effects that we ascribe to disclination defects J18j . 

• The 2"*^ step is to exploit the absence of mass in the action ([3]), that, in the case in point of 
the Dirac field, guarantees full local Weyl symmetry [3]. In the case of conformally flat metrics, 
gfiu = (p'^Vfiu, Weyl symmetry gives 

i f d^^x^i^-f'^Vf^iP = i f ct^xi^'-i'^da^' (4) 



where, under Weyl's transformations, ip = (j)~^ip' . Namely, the flat space theory is classically 
equal to a curved space theory. We have to consider quantum effects, though. As we are in a 
odd-dimensional spacetime, this does not refer to quantum trace anomalies, absent here |19j . 
What we refer to here (see [ll|20]) is that one can, essentially, do two type of things: (a) perform 
measurements in a frame where the quantum effects of curvature disappear, e.g., 

S{xi,X2) = {0\iP(xMx2m = '{0\iP'{x^)i;'{x2)\0y = S'{XI,X2) , (5) 

where S{xi,X2) is the two-point function in the curved case, and S'{xi,X2) is the two-point 
function in the flat case; or (b) perform the measurements in a frame where, due to Weyl 
symmetry, the effects are seen as a simple modification of the flat case (conformal triviality [5j ) 

S{xi,X2) = '(0|V'(xi)Vi(x2)|0)' = (^-\xi)<p-Hx2) S'{xi,X2) . (6) 

In our case, due to the obvious difficulty to realize a frame where the measuring apparatus is in 
a gravitating field whose spacetime coincides with the graphene curved spacetime, the second 
choice will be the right one. 



• The 3"^ step is to stick to a metric that is the most easy to reahze in a lab, namely 

(2) 

i.e., with all the curvature in the spatial part g^^a, and no time-dependence. Can this be 
conformally flat, rich and still interesting? To answer all these questions we need to solve 

PIE] 

Since all surfaces of constant Gaussian curvature /C solve this equation |lj, they give rise to 
conformally flat (2-|-l)-dimensional spacetimes g^iu ■ Hence the answer to the above questions 
is yes: by simply curving the graphene sheet one can have rich enough situations (surfaces of 
constant /C) implying a conformally flat spacetime. 

• The 4**^ step. The latter result is intrinsic, but to implement Weyl symmetry we need the 
explicit conformal factor, and this might mean to change the frame 

(?'^ = (t,x,y)^Q^ = (r,X,y) (9) 

so that 

g%^'\Q) = <i>\Q)9lf{Q)- (10) 

This can be technically difficult (a system of 6 partial differential equations), but the most 
important here are the physical meaning, global predictability, and practical feasibility of this 
coordinate frame. 

The surfaces of /C = -|-l/r^= const, can be described by the spatial line element |21j 



d£2 ^ ^^2 ^ ^2 cos2(u/r)du2 (11) 

where u is meridian coordinate, and v the parallel coordinate, and with the real constant c such 
that c = r (sphere), c < r^ c> r. These Q^ are difficult to find. 

Interesting objects, like, e.g., the black-hole of [H], are related to spacetimes of negative 
constant curvature. For the surfaces of /C = — l/r^= const, the following theorem by Hilbert 
holds: No non-singular surface of constant /C = — 1/r^ can he embedded in R^. If we can 
live (and we can) with singularities, such surfaces are many more than the positive curvature 
counterparts (there are, actually, an infinite number of them). Among them there are the surfaces 
of revolution, called pseudospheres, and they can be described by the spatial line element [21] 



di^ = du^ + (c\ cosh^{u/r) + cl sinh^ (u/r)) dv"^ (12) 

with Cl, C2 G R. For all these surfaces (and also for those not included in the above line element) 



"^graphene "'^ r2l"3; -j- Ciy ) -^ 



—dt — dx — dy 



(13) 



with X, y the abstract coordinates of upper-half plane model of the Lobachevsky geometry (see 
[16] for the case in point of graphene, and [22] for a general introduction to the subject). To make 
sense in a laboratory, these coordinates need to be expressed in terms of measurable coordinates, 
say the (n, v) coordinates above, and this amounts to choose the surface. The line element in 
square bracket is flat, but this is not yet the solution of our problem, because we need to see, 
case by case, to which physical Q^ the coordinates (t, x, y) correspond to. 




Figure 4. The R^ coordinates of the Beltrami pseudosphere, in the canonical form, are [23] 

x{u,v) = R{u) COS V, y{u,v) = R{u) sinv, z{u) 

with R{u) = ce"/^ c > and r = V-ZC" 

curvature. We can choose c = r, thus R{u) £ [0,r] as u G [— oo,0]. The surface is not defined 



r(y/l — R^{u)/r'^ — arctanh-yi" 
> where /C is the constant negative Gaussian 



R\u)/r^ 



for R > r {z becomes imaginary), an instance of the Hilbert theorem. In the plot r 
uG [-3.37,0], V e [0,27r]. 



1 and 



As a counter-example, consider the case of the hyperbolic pseudosphere (ci 7^ and C2 = 
in ([ISD) 

X = e'^^/'^ tanh(u/r) , y = e^^V cosh(n/r) . (14) 

Since the conformal factor r"^ /y"^ in (J13p is multivalued {v is an angle), this solution as no global 
predictability in q^^ = {t,u,v), hence one must look for a different frame where this does not 
happen, but the frame in point might be difficult to realize in practise (and surely it is not the 
lab frame q^, where all that needs to be done is to curve the surface). 
It is the Beltrami pseudosphere (see Fig. H]) that solves the problem 



v/r , y = e ^'^ /r 



(15) 



i.e. r jy makes sense over the whole surface/spacetime. Hence, the physical Q^ are the lab's 
q^ . A physically doable 5^57'''^^'^*^ = 0^ 9^ is already there in the frame q^: 



ds: 



;raphene 



fiu 



ds 



(B) 



^2u/t 



-^-/^■{dt'-du^) 



r dv 



(16) 



Now we can say that the line element in square brackets is the Rindler line element, ds?j^s, and 
exploit this fact in our quest for measurable effects of the QFT in curved space interpretation 
of the physics of graphene vr electrons. 

• The 5*^ step consists in implementing the Unruh requests on this setting. In particular it 
is important the choice of the quantum vacuum of reference, taken here to be the Minkowski 
vacuum |0)' = \0)m, thus the Green's function of interest are of the form (see ^ and [HEjfTU]) 
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S"{xi,X2) = {^M\^"{xi)^"{x2)\QM)=<i)'\xi)(t)~'{x2) {^M\i^''{xiW{x2)\^M) 



(17) 



p(1.5xl0i^eV-^m-2) 

l.Qr 



0.8 



p<'*'(u=-r) 



p<*'(u=-3r) 



r= 10;um^T'o~0.13K 
O.d 0(^T'~0.13K)<u<-30^m(^r-0.05mK) 




-0.10 



-0.05 



0.00 



0.05 



0.10 



E (meV) 



Figure 5. Series of dl/dV spectra ("line-cut") from an STM. The model applies to energies 
below the natural scale, E* = Uvf/t — 6.6 x 10"'^ eV x [r mm]""'^, where [r mm] is the numerical 
value of r measured in mm. In this plot E E [— 0.1,0. l]meV, corresponding to r = 10/um. The 
flat LDOS is shown in blue, u g]0, — 3r], u = (red dashed curve) corresponds to the singular 
boundary of the pseudosphere, beyond the upper limit u = — 3r (green curve) the pseudosphere 
is too sharp. To each spectrum corresponds a u-dependent temperature T while the constant 
temperature, same for the whole line-cut, is To = e~'"'^T ~ 0.13K. 

where B stands for "Beltrami", R for "Rindler", and the Green's function on the right side is 
S^, the customary two-point function used in the computation of the Unruh effect. The power 
spectrum is customarily computed by taking S^ at the same point in space, q, and at two 
different times, with t2 — ti = e^'^'r, and r the proper time (the factor e"'^ is related here to the 
proper acceleration, and needs be constant, to fit the requests for the standard Unruh effect). 
From this, by Fourier transforming on the time variable, we get [15] 



?B.^ 



I. 



F"{u,q)^-T, 
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dre-^-^5«(r,g,C) 



(18) 



and it coincides here with the (not yet physical) electronic local density of states (LDOS) [9] 



vr 
As this is a massless case, the result of the computation is exact 
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(19) 



(20) 



where T = a(u)/(27r) = e^'^ /{2iTr). We are now in the position to make a verifiable prediction, 
by going back to the physical spacetime. Due to the Weyl symmetry, this is easy (see (|17p ) 



Ri 



p^(^,gl = 0-^(g>«(t^. 



(21) 



Thus, the LDOS predicted for a graphene sheet shaped as a Beltrami pseudosphere is [2] 



^ ,41 Se-2"/'- 



P ^^' ^' "■) - vr {UvF? exp [i?/(A:Broe-A)] _ i ' (^2) 

i.e., a Hawking phenomenon with Hawking temperature 

ro-:^f^. (23) 

Zti Kb r 

In Fig. [5] a plot of what we expect from a dedicated experiment with a Scanning Tunneling 
Microscope [2]. 

3. Paving the way to future experimental set-ups 

Lately, I have been interacting with experimentalists from the Manoharan Lab of Stanford 
University [23], who asked whether it is possible to see effects, of the kind described above, 
on their "molecular graphene" [25j, that is a planar hexagonal lattice of carbon monoxide, 
CO, molecules that mimics graphene. Unfortunately, there were severe limitations in their 
experimental settings. The first is the fiat setting, the second is the small number of points, the 
third is that molecular graphene is a simulator of the real thing that, although a very good one 
to test certain properties, it has basic parameters (e.g., the neighbor-to- neighbor distance, i) 
that make our task more difficult. 

To write things in a list, what I had to address are the following issues: 

(i) How to tile a Beltrami pseudosphere? 

(ii) Would a "planar development" (the result of mapping the surface on a plane, hence with 
curvature traded for strain), still captures the features necessary for the Hawking-Unruh 
phenomenon to happen? Which one, among the infinitely many, is the best? 

(iii) Is the effect still there for a generic surface of negative constant /C? If yes ... 

(iv) ... can we solve issues (i) and (ii) for the given surface? 

(v) Can we spot any Hawking-Unruh in a planar configuration of molecular graphene, that 
mimics the planar development of a portion of a generic graphene surface of (not even 
necessarily constant) negative /C? 

My conclusions for the last point (the best meeting point between Manoharan Lab's and this 
model's needs), are negative, as I will explain below. This analysis gave extra arguments in favor 
of a dedicated experiment, with a Beltrami pseudosphere, as the best way to test this model. 
Nonetheless, the attempt was worth the effort. First, I learned that point (iii) has a positive 
solution. Thus, the realization of a generic graphene surface of constant negative Gaussian 
curvature, although not being the best possible setting, should be enough to test a modified 
formula (j22p . at least in an approximate way. On this I will briefly comment below, and will 
write in detail elsewhere. Second, for whatever dedicated experiment one would like to carry 
on, matters like the issue (i) in the list are crucial. In the following, I will address in detail this 
latter issue, and will also briefly comment about the other issues. 

Therefore, let me focus on the geometric and topological constraints on building a Beltrami 
pseudosphere with the hexagonal graphene membrane. There are two issues here. The flrst is 
the inflnite (non-compact) surface (u G [— oo,0]). The second is how many heptagonal defects 
are necessary. Let me introduce, briefly, the matter by recalling how to tile a compact surface 
S (see, e.g., [26j). Suppose we want to tile S with regular s-sided polygons assembled in such a 
way that each vertex is shared by 3 polygons, and each edge is shared by 2 polygons. If Ug is 
the number of polygons used (there are, say, n^ pentagons, ng hexagons, n^ heptagons, etc.) to 



tile S, the resulting polyhedron P has Vp = ^/"i^n^sS vertices, Ep = l/2^^nsS edges, and 
Fp = ^„ Ug faces, giving for the Euler characteristic 

X(S) = Vp-Ep + Fp, (24) 

the following expression 

5]n,(6-s) = 6x(S). (25) 

n 

For the sphere, the Gauss-Bonnet theorem [,23J, that links the topological number, Xi to the 
geometric object that is total curvature, gives 

Ktot = JlCiffi = 2TTx, (26) 

and, being /C = 1/r^ = constant 



X = 7^-3' (27) 



Area 
27rr2 

or, being Area = 47rr^, we have x = 2. Hence, if we try to wrap the hexagonal lattice of graphene, 
to form a tiling of the sphere, we cannot match J2n ^s(6 — s) = 12, only with hexagons (whose 
number is not constrained by this argument), but need something else. Considering that the 
creation of each defect costs energy, the minimum number of defects we need are 12 pentagons, 
X]„"'s(6 — s) = 77.5 (6 — 5) + uq (6 — 6) = 12, with which we could make the smallest (and 
"minimum energy") tiling of the sphere, the dodecahedron (F = 12, V = 20, E = 30). Adding 
20 hexagons to it, one makes a buckyball (fullerene). One can increase ng further (according 
to the rule ne = 10(r - 1), with T = m^ + P + ml = 0,1,3,4,7,9,..., and excluding T = 0), 
making "inflated" buckyballs. Once we trade the continuous sphere for its tiling polyhedron, 
the curvature resides in the vertices. For the sake of simplicity, let us focus on the icosahedron, 
dual {F 'ir^ V, E the same) to the dodecahedron, so that Ktot = 12 Ky" = A-n, and Ky° = K^. 
The curvature associated to each pentagon is then 

i^5 = +^/3 , (28) 

while no curvature is associated to the hexagons, Kq = 0. Of course, the same results can be 
obtained using the dodecahedron, although the argument becomes more involved. A direct way 
to obtain this result is with the defect/excess angle argument. The curvature at one vertex of 
the polyhedron is 

Kv = 2T:-Y^ai, (29) 

i 

where aj is the angle between edges of each face at the vertex. For a pentagon (formed with 
equilateral triangles in the icosahedron), ai = 7r/3, and i = 1, .-.,5, hence, Ky = 211 — 57r/3 = 
-|-7r/3 (defect angle). The same argument, for an hexagon gives Ky = 2tt — Qtt/3 = 0. 

Notice that the radius of curvature, r, does not enter these formulae, hence the curvature 
associated to each defect is fixed. If we want to increase the curvature of the sphere (make r 
smaller), we do not have to increase the number of defects, but their density (i.e., we have to 
decrease the number of hexagons). 

From the elastic point of view, it has been shown that, for the sphere, the creation of 
heptagonal defects is energetically favored in certain conditions ^27j , but this does not violate the 
Euler- Poincare result (^^, provided they come in pairs with pentagonal defects, n^ — n^ = 12, 
i.e. n^ = 12 + X and n-j = x^ and that is indeed the case (scars [27]). The same argument 
descending from the Gauss-Bonnet theorem (j26|) (and the excess angle argument of (j29|) ) gives 

Kj = -7r/3 , (30) 



as the curvature associated to each one of these defects. Thus, when we want to build (tihng 
of) surfaces of negative K. we need heptagonal defects. 

If we want to wrap the graphene membrane into a Beltrami pseudosphere, we first need to 
consider that the issue of an infinite (non-compact) surface is actually minor, from a practical 
(lab) point of view. One can truncate the Beltrami, without changing the topology, and still 
has a surface of constant negative curvature. (Notice, also, that the elliptic and hyperbolic 
pseudospheres are finite surfaces, so, for them, this problem is absent). So one can think of 
tiling the surface in the standard way illustrated above. The issue here is the number and types 
of defects necessary, and, in turn, how the Euler-Poincare and Gauss-Bonnet theorems apply 
here, as they have been conceived for compact surfaces. 

Having the topology of an annulus, the Euler number for Beltrami is 



XBeltrami = . (31) 

(The same is true for the other two pseudospheres, not clear is the matter for more exotic 
pseudospheres [28]). If one blindly applies (pl|) . V — E + F = — )• X](6 — s)ns = 0, the 
implications do not make sense. Would that mean that s = 6? So I can tile with hexagons 
all over? But that is true for the plane, the cylinder, or the cone, not for the pseudosphere. 
Would that mean that 71,7 = ns? But then I can satisfy the constraint with even 1 heptagon 
and 1 pentagon, and this makes no sense either. The unsolved (and intriguing) problem here is 
that, as far as I know, it is not known the negative curvature counterpart of the dodecahedron 
(or, of its dual, the icosahedron), the "elementary" polyhedron one obtains from the Euler- 
Poincare for the sphere. There must be a fixed number of heptagons, that take into account 
the basic geometric properties, just like the 12 pentagons do for the sphere. The beauty and 
intricacies of the Lobachevsky geometry come in here, and make changes to the sphere story, and 
then introduce the singular boundaries, but the basic object should be there to give a unifying 
"Platonic-Lobachevsky" solid for all the infinite pseudospheres (not just the three that clearly 
have X = 0) but all the others). A partial solution to this problem, and in turn to our main 
problem (wrapping graphene to make a Beltrami), is to consider that the Gauss-Bonnet for 
surfaces with boundaries is [23] 

/ /CdV + I ICgdl = 2ttx , (32) 

JE JdT. 

with K.g the geodesic curvature that measures how far is a curve from a geodesic (e.g., /Cg = for 
the equator of a sphere). In (|32p there are two parts, one for the bulk that is the total curvature, 
Ktoti and one for the boundaries. It is Ktot that tells about the number of defects, the boundary 
term does not. Indeed, from Area = 27rr^, we have Ktot = — Area/r^ = — 27r, (the same is true 
for the other two pseudospheres), hence, using ([301) . Ktot = — 27r = ny K-j = n-j (— vr/S), that 
gives 

n7 = 6 . (33) 

In other words, the Euler-Poincare and the Gauss-Bonnet theorems "decouple" for 
pseudospheres, and we cannot apply (p^ to compute the number of defects, but need the 
curvature arguments just given. That said, though, everything else is like for the sphere: the 
negative curvature counterpart of the icosahedron/dodecahedron, at least for the three classic 
pseudospheres, has 6 heptagons, arranged to form a non-compact "polyhedron" (actually, I did 
not succeeded to use only heptagons, and a certain number of hexagons are necessary even for 
the basic object), and, elastic energy allowing, one can add pentagon- heptagon pairs. This is 
my solution of the issue (i) in the list above. Let me now face the other issues. 

The procedure to map a surface onto a plane is called "planar development" ^3]. For the 
intrinsically fiat surfaces, this can be done isometrically (that is why those surfaces are called 



"developable"). For intrinsically curved surfaces, planar development comes at the price of 
deforming the distance among points (think of any map of the Earth), hence of trading curvature 
for strain, and there are infinitely many different "planar versions" of the curved surface, see, e.g., 
[29]. At first sight, the effects we are describing with (|22p. being due to intrinsic curvature, could 
never be seen on a planar development. In the case in point of the graphene membrane, though, 
there is some hope to still have some sort of signal of the Hawking-Unruh effect, although we can 
expect it to be substantially deformed by the planar development. It is so, because, in our model, 
the intrinsic curvature of the graphene membrane is encoded in the heptagonal defects, hence, 
retaining these defects in the developed Beltrami surface, in principle, could be the way out to 
test the effect on a flat sample. Unfortunately, the geometric constraints on the development, 
and the limits of the experimental set-up of the device of Manoharan Lab, conspired to make 
this test impossible. Let me explain in a bit more detail. 

The actual values of various important quantities differ, for graphene and for molecular 
graphene (for the latter I rely on [31] , see also [24J and [25] ) . In the table [H I summarize the 
ones relevant for this discussion. There, I give two energy scales, E^ ~ fivp/t (below which the 
continuum Dirac field approximation holds), and E^ — fivpjr (below which the field theory in 
curved space time holds). The formula (|22p has a chance to be tested only if r >> £, and it 
holds for energies below E^. On this see [3] and [2j|T6]. Essentially, what goes wrong with the 
Manoharan Lab set-up is that, even to just include all the necessary 6 heptagons, the limited 
number of points makes the density of defects (hence, the corresponding curvature) too close to 
the limit for the model to hold. Given the values of vp and £, and the maximum sensitivity of 
the STM machine of 1-10 meV, to a small enough curvature (big r) corresponds an energy too 
small to detect the effect in that range (see the entry for r ~ 100£ in table [U and consider that 
the electrons need to have an energy below the corresponding E.^). 



Table 1. Values of various important quantities, for graphene and for molecular graphene. Here 

Vp = ^t^, with t hopping parameter, i nearest-neighbor distance, and Ei ~ %vp/l, E^ — hvp/r. 





Graphene 


Molecular Graphene 


i{m) 


1.4 X lO-^'^ 


20 X lO"^'^ 


vp (m/s) 


9x 10^ 


6.5 X 10^ 


Ei (eV) 


4.2 


0.2 


Erir ~ 100£) (meV) 


42 


2.1 


Er{r^50i) (meV) 


84 


4 


Er{r ~ 20i) (meV) 


200 


10 



Indeed, the number of defects per unit area is S"^' = n-jj A = S/vr r^^, as 77.7 = 6 and A = 27rr^ 
for Beltrami. So the number of defects in the area a of the surface is n" = dS ' a = S/vr a/r"^, 
hence 



r{n'^) 



2>a 



TTUy 



(34) 



The number of CO molecules is limited to 500, hence an estimate of the area is: 170 hexagons, 
each with area 3\/3/2£2^ giving a ~ 450^^ ~ 9000A^. With these, Er ~ hvp/r = a~^E^ = 
a~^ X 0.2eV, where a = r/£, all the quantities are for molecular graphene. In table [21 I report 
the estimate of the values of r as a function of the defects present in the lattice, and the 
corresponding E^s. It is clear from there that we are right at the limit of validity of the model, 
(see the entries for Uj = 6). 

Furthermore, a planar development of the Beltrami, without any meridian cut, is a highly 
complicated problem, as is well known to shoes, ships and aircrafts makers [32]. Our goal was 



to "squash" it into a deformed annulus [30], without cutting the surface along a meridian (the 
truncation along some parallel on the "tail", instead, was required, as explained earlier). This 
would help retaining the necessary global properties of the Beltrami spacetime. Unfortunately, 
we realized that the strain would be too big [30j, and would break the cr-bonds [3lj . 



Table 2. Estimates of radii of curvature as function of the number of defects, and corresponding 
threshold energies, for a 500 points lattice of molecular graphene. The last values correspond to 
include all the six heptagonal defects, required to mimic a Beltrami (and other) pseudosphere. 
The corresponding curvature is too close to the limit to have reliable results. 



n? 


rji 


Er (meV) 


1 


21 


10 


2 


15 


14 


3 


12 


17 


4 


10.5 


20 


5 


9.5 


22 


6 


8.5 


25 



For all these reasons, the last proposal that was left to test the model, within these limits, 
was to see whether some relics of the effect are still there, for a configuration of planar molecular 
graphene that includes few heptagonal defects (with this, at least the condition of small curvature 
of the original surface is better satisfied/less unsatisfied, see table [2]). In this fashion, the sample 
could mimic the planar development of a 'portion of a generic (not necessarily constant) negatively 
curved surface. The answer, as said, is negative. Let me explain why. 

One piece of good news here arrived from the discovery that the issue (iii) of the list has 
positive answer. Namely, it is indeed possible to have a modified and approximate version of 
formula ()22p at work, for small neighbors of a point of a graphene membrane shaped into a 
generic surface of constant negative curvature. Indeed, it is true that, as discussed in the 5*^ 
step of the previous Section (see [2] and [T6] ) , that the global predictability for the given surface 
in the lab frame is lost, but, as long as local considerations are concerned, some predictive power 
of ()22p is still at work. Full details of this result will be given elsewhere. For the moment let 
me say that the main reason for this is that the horizon, in this model, is reached "time-wise", 
like in the standard Unruh effect, and not "space-wise", i.e. the horizon is not a specific curve 
on the manifold. Being the time t the same for all surfaces in the lab frame q^ ^ the effect is still 
there for all, although modified and only valid portion by portion of the surface. 

Nonetheless (leaving aside the issue of constant vs nonconstant /C), one thing is to focus 
on one portion (neighbor of a point) of a whole surface, for which the global conditions for the 
existence of an horizon (hence the Hawking) are realized, and another thing is to have a detached 
(cut-away) portion of that surface. For instance, in order to feel the curvature, the electrons 
need ho have enough time to travel back and forth on the whole surface, an instance that also 
means that the electrons that feel the curvature are those of small energy. When this takes place 
on a surface, the electrons come back to the initial point (for a non-compact surface, at least 
some of them), while for a piece cut-away from the surface, obviously this cannot happen. Thus, 
there is no hope to see any effect that originates from global properties, as the Hawking-Unruh 
of (I22|). 

To conclude on a positive note: although it was impossible to test these ideas in the 
experimental settings of Manoharan's lab (that were not meant to test these kind of effects, 
in the first place), the fact that we can test them on a generic surface of constant negative /C is 
very promising for future, dedicated experiments. 
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